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1.0  INTRODUCTION 

ANCO  Is  investigating  the  use  of  radiographic 
techniques  to  deteraine  properties  of  anatoalcal  segaents  - 
aass.  center  of  gravity,  and  the  inertial  tensor.  As  shown 
in  our  proposal,  three  orthogonal  projections  (x-ray 
photographs)  through  the  segnent  will  deteraine  these 
properties  if  sufficiently  high  energy  x-rays  ere  used  so 
that  attenuation  ia  independent  of  atonic  nuaber  and  Juat 
dependent  on  aass  density.  This  occurs  above  a  few  hundred 
keV  energy. 

The  intent  of  this  study  is  to  investigate  the 
practicality  and  potential  accuracy  of  the  aethod.  Progress 
to. 'gate  is  described  herein.  All  work  will  be  perforned  by 
the  end  of  February  and  the  work  plan  is  described  here  as 

I 

2.0  X-RAY  SOURCE 

The  decision  aust  be  aade  as  to  the  source  of  the  / 
x-rays,  either  an  x-ray  tube  can  be  used  or  an  x-ray  (gaaaa 
ray)  source.  Most  conaercial  x-ray  tubes  operate  below 
several  hundred  keV.  Special  industrial  x-ray  tubes  can 
reach  several  MeV.  All  tubes  provide  a  distributed  energy 
source  containing  energies  above  and  below  their  specified 
energy  (the  bulk  being  below). 

X-ray  (vs.  gaaaa  ray)  sources  can  provide  single  energy 
x-rays  at  high  energies.  Exaaples  Include  Ceslua  137  at  .60 
NeV  (33  year  half  life.  .39  r/hr/Curie  at  P  aeter  dosage) 
and  cobalt  60  (9.3  year  half  life.  1.35  r/hr/Curle). 
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Our  studies  indicate  that  a  radioactive  source  is 
preferable  to  an  x-ray  tube  for  the  following  reasons: 

•  The  source  can  have  higher  energy  than  the  nore 
coaaonly  available  tubes. 

•  The  source  energy  is  aonocroaatic  (single  energy) 
while  the  tube  has  auch  radiation  at  lower  energies 
where  absorption  is  atonic  nuaber  dependent.  This 
is  true  even  for  high  energy  tubes. 

•  Sources  are  less  expensive  and  saaller  than  tubes, 
asking  the  exposure  set  up  less  expensive  and 
slapler.  Also  the  saaller  sources  can  be  placed 
closer  to  the  subject  without  blurring  due  to  source 
size  and  are  easier  to  shield. 

•  The  source  is  nore  easily  shielded  than  is  a  aachine 

because  of  the  snail  size  of  the  source. 

3.0  EXPOSURE  TINE  AND  DOSE 

Prellalnary  calculations  indicate,  for  exaaple,  for 
Agfa  D-7  fila  (with  5  ag  Pb  screens  and  X10  lntenaifler)  and 
the  source  placed  1  aeter  froa  the  segaent  and  screen,  the 
following  exposure  tine  and  doses  would  be  required  (the 
source  strength  below  one  of  those  currently  available  to  us 
and  could,  of  course,  be  varied  in  future  work). 


Source 

Strength 

Exposure  Tiae  (ain) 

Dose  (roentgens) 

Csl37 

.146  Ci 

343 

2.2 

Co60 

.80  Cl 

26 

0.6 

The  exposure  tine  can  be  reduced  by  placing  the  source 
closer  to  the  object  or  increasing  source  strength  (dose 
will  not  be  significantly  changed).  Thus,  for  exaaple.  use 
of  a  one  curie  source  at  1/2  aeter  would  result  in  an 
exposure  tiae  of  12.5  ainutes  for  Csl37  and  5.2  ainutes 
for  Co60.  Both  this  distance  and  source  strength  are 
practical.  As  several  exposures  (typically  3)  aust  be  Bade. 


we  are  looking  at  total  exposure  times  of  15  -  40  minutes,  it 
appears  feasible  to  hold  an  anesthetized  aniaal  still  for 
this  length  of  tlae.  A  dose  of  2-7  roentgens  should  not 
cease  the  enleel  here  on  a  one  or  two  tine  basis.  It  is 
possible  that  alternate  files  eay  further  reduce  these  tlaes 
and  dosages.  Larger  sources  could  be  used  to  reduce 
exposure  tiee. 

4.0  EXPOSURE  CONFIGURATION 

/ 

As  anticipated,  eultlple  exposures  are  desired.  The 
Issues  to  be  considered  are: 

•  Can  non-orthogonal  views  be  used? 

•  Can  aore  than  3  views  be  useful? 

•  Can  exposure  be  aade  simultaneously  using  eultlple 
sources? 

•  How  can  the  anatoalcal  segment  be  "spotted"  to 
establish  a  coordinate  reference? 

The  first  three  points  above  have  a  coenon  solution. 
Non-orthogonal  views  are  useful  so  as  to  better  Isolate  a 
given  segment.  For  example,  the  head  is  connected  to  the 
body  by  the  neck.  If  the  aass  properties  are  desired  above 
a  given  reference  (say  a  specified  vertibrae)  and  if  one 
could  severe  the  head  at  that  point  and  "float"  it  in  space, 
one  could  easily  take  3  orthogonal  views,  each  encoapassing 
the  entire  head.  Since  we  wish  to  do  non-destructive 
testing,  this  can  not  be  done.  Orthogonal  views  of  an 
unsevered  head  will  result  in  soae  non-coverage  and  error, 
as  illustrated  in  Figure  1. 

Non-orthogonal  views,  as  shown  in  Pigure  1.  reduce  this 
error.  The  error  could  be  reduced  to  zero  by  having  the 
views  in  the  sue  horizontal  plane.  However  they  would  no 
longer  be  independent  and  this  would  lead  to  an  inability  to 
estiaate  certain  aass  properties.  (Reaeaber  that,  as  shown 
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in  our  proposal,  two  orthogonal  views  in  a  horizontal  plane 
could  deteraine  aass.  center  of  gravity,  and  all  but  the 
off-diagonal  teraa  of  the  Inertia  tensor .  Reaeaber  also 
that  If  aany  views ,  each  differing  by  a  saall  angle,  were 
taken  la  a  horizontal  plana,  we  would  have  a  "CAT"  scan  and 
could  deteraine  the  density  distribution  and  hence  all 
properties.  He  are  however  trying  to  deteraine  all 
properties  with  a  few  scans.)  Hence  soae  angle,  however 
saall,  aust  exist  between  all  scans.  Multiple  scans  produce 
redundancy  and  help  to  aake  up  for  alaost  dependent  scans 
and  are  consequently  potentially  useful. 

In  addition  to  these  considerations,  we  oust  concern 
ourselves  with  beaa  spreading  (parallax,  fan  beaa)  as  we  plan 
to  use  a  source  close  to  the  segaent.  A  source  closer  than 
about  2  aetera  will  have  a  significant  non-parallel  beaa. 
Non-orthogonal  views,  aultiple  views  (aore  than  3)  and  beaa 
speed  all  introduce  transforaations  of  the  data  that  auat  be 
accounted  for.  Fortunately  there  is  a  staple  technique  for 
doing  so. 

Assuae  that  any  nuaber  of  views  are  taken  in  any  nuaber 
of  orientations.  The  views  are  discretely  converted  to  give 
a  vector  R  of  aeasureaents .  For  exaaple,  if  5  views  were 
taken  and  a  20x20  ca  grid  used  for  each,  then  there  would  be 
5x20x20  or  2000  aeasureaents,  and  R  would  be  a  vector  of 
length  2000.  Assuae  also  that  in  a  convenient  orthogonal 
coordinate  systea,  the  space  containing  the  segaent  is 
discretely  converted  into  cells  of  assuaed  unifora  density 
P.  Thus  if  the  space  is  20x20x20  ca,  we  can  define  a  vector 
P  of  size  8000.  Depending  on  the  geoaetry  of  the  views  and 
the  spreading  ray  path,  each  eleaent  of  R  is  linearly 
dependent  on  soae  (in  fact  a  very  few)  of  the  eleaents  of 


where  T  is  a  large  (2000  x  8000)  matrix  with  most  terms 
equal  to  zero.  The  P  can  not  be  solved  for  as  there  are  too 
few  equations  (2000)  for  the  unknowns  (8000)  and.  further, 
all  the  equations  may  not  represent  linearly  independent 
data.  We  do  not.  however,  wish  to  know  P  but  rather  certain 
linear  coebinationa  of  P  (l.e.  aass.  center  of  gravity, 
inertia  tensor  -  see  equations  1-4  of  our  proposal).  For 
exaaple  the  elenents  of  the  inertia  tensor  can  be  placed  in 
a  six  coaponent  vector  J  and  equation  3  and  4  of  our 
proposal  written  as 

J  -  UP 

where,  for  exaaple  U  is  a  6  x  8000  aatrix.  For  the  noaent. 
assuae  that  the  aatrix  T  can  be  inverted  and  call  this 
Inversion  TI..  Then 

P  -  TIR 

and 

J  -  UP 

or 

J  -  UTIR 

where  UTI  is  a  6  x  2000  aatrix.  We  would  in  fact  have  an 
over  constrained  systea  (aore  equations  than  unknowns). 

Note  that  UTI.  if  it  existed,  is  only  a  property  of  the 
geoaetry  of  the  sources  and  views  and  is  independent  of  the 
properties  of  the  segnent  under  study. 

This  fora  of  problea  occurs  frequently  in  various 
processing  and  estiaatlon  tasks.  In  fact  a  aatrix  TI  can 
be  defined  which  has  aany  of  the  properties  of  the 
(non-existent)  inverse  of  T.  It  is  called  the  Pseudo 
Inverse  (also  the  Penrose  inverse,  or  singular  value 
decoapositlon  inverse).  Without  further  elaboration  here. 
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note  that  J  =  L'TIR  should  provide  a  good  estimate  of  the 
inertial  properties.  We  know  that  in  an  orthogonal  parallel 
scan,  the  answer  is  exact  (aa  provided  in  our  proposal)  and 
feel  that  non-orthogonal  scans  with  spreading  beans  should 
also  yield  a  good  estlaate  (aa  long  as  the  angles  are  not 
too  far  froa  90  degrees)  and  that  aultlple  views  (beyond  3) 
can  help.  Thus  the  pseudo  Inverse  allows  solution  In  one 
foraalisa,  accounting  for 

•  Non-orthogonal  views 

•  Spreading  beans 

•  Redundant  views  (in  the  Halt  -  a  CAT 

Scan) 

The  evaluation  of  UTI,  while  a  lengthy  conputatlon. 
needs  only  be  done  once  and  will  not  change  once  a  view 
geoaetry  is  chosen. 

These  are  two  other  Issues  nentloned  above  - 
slaultaneoua  exposure  and  spotting  for  coordinate 
reference.  In  principal  all  views  could  be  exposed  at  once, 
thus  reducing  the  tlae  the  anlaal  aust  be  held  still 
(anesthetized).  There  is  no  aatheaatlcal  reason  why  the  data 
could  not  be  analyzed  even  if  one  view  received  direct 
exposure  froa  aore  than  one  source.  However,  the  effect  of 
scattering  would  be  greater  in  such  a  situation.  This  is 
because  with  one  source,  nost  of  the  scattered  radiation  is 
directed  away  froa  the  fila.  With  aultlple  slaultaneoua 
sources  and  views,  aost  of  the  scattered  radiation  would 
fall  upon  one  of  the  views.  This  effect  will  be  quantified 
but  will  probably  suggest  against  siaulataneous  exposure. 

Note  that  if  used,  slaultaneous  exposures  reduce  exposure 
tlae  but  not  dosage. 

The  spotting  of  a  coordinate  reference  is  required. 

Our  idea  is  to  place  a  snail  dense  object  (ball  bearings)  at 
known  points  on  the  segnent  (e.g.  vertebrae,  top  of  head. 


temples,  etc.)  and  then  locate  these  in  each  of  the  view 
exposures.  The  data  reduction  computer  program  would  then 
relate  these  to  the  coordinate  system  in  which  the  mass 
properties  are  defined.  The  weight  of  the  objects  would 
have  little  effect  on  the  calculated  mass  properties. 

S.O  ACCURACY 

A  primary  goal  is  to  detereine  the  potential  accuracy 
of  this  technique.  Note  that  researchers  using  a  CAT 
scanner  have,  at  great  cost,  produced  results  with  a  5* 
accuracy.  It  Is  unlikely  that  the  accuracy  of  our  few  scan 
techniques  would  be  better. 

There  are  many  factors  affecting  accuracy. 

•  Discretization  -  A  coarser  grid  will  give  greater 
error  than  a  finer  grid.  Ve  are  tentatively  working 
with  a  1  ce  grid.  Coaputer  sinulatlonn  will  be  used 
to  suggest  the  relative  errors  between  a  1  ce  and  a 
cm  grid. 

•  Attenuation  sensitivity  -  Even  at  several  hundred 
keV  and  higher,  there  is  soae  variation  in 
attenuation  that  depends  on  atonic  number  (not  Just 
density).  Computer  simulations  will  also  be  used  to 
suggest  the  effect  of  this  error  for  typical 
variations. 

•  Scattering  -  Not  all  the  radiation  is  either  passed 
or  absorbed.  Some  is  scattered.  If  the  scattered 
radiation  did  not  reach  the  film  plate,  it  could  be 
combined  with  the  absorbed  radiation  for  an 
effective  absorption.  However,  soae  of  the 
scattered  radiation  nay  reach  the  film,  perhaps 
after  multiple  scattering  and  perhaps  downgraded  in 
energy.  Our  initial  computations  suggest  that  about 
1M  of  the  radiation  reaching  the  film  nay  be 
scattered  radiation.  Moving  the  film  away  from  the 
source  and  segment  can  reduce  thia  fraction.  It  is 
not  yet  clear  how  much  error  this  10*  will  introduce 
but  it  is  felt  to  be  less  than  10*.  The  experiments 
planned  should  suggest  the  magnitude  of  the  error. 

•  Calibration  and  Dynamic  Range  -  All  films  have 
certain  sensitivity  and  dynamic  range.  The 
sensitivities  are  non-linear  functions  of  exposure 


and  have  a  threshold  and  saturation.  Thus  too  snail 
an  exposure  will  produce  a  zero  signal  regardless  of 
variation  of  the  density.  Too  large  an  exposure 
will  produce  a  "100*"  signal  regardless  of  the 
variations  of  density.  A  fila  and  exposure  tine 
aust  be  chosen  to  produce  a  result  in  the 
quantitative  region  of  sensitivity  and  the  resulting 
radiograph  nuat  be  capable  of  being  quantitatively 
read. 

Initial  slaulatlona  suggest  that  the  attenuation  of 
the  bean  will  vary  between  0*  and  90*  (leaving  100* 
to  10*  of  the  bean)  .  Thus  a  dynaalc  range  of  about 
10  la  required.  The  aaxinun  attenuation  is  expected 
to  equal  about  1"  of  steel  (the  approxlnate 
equivalent  path  density  of  a  ray  through  the  head). 
To  provide  calibration  of  the  file,  a  1"  steel  and 
1"  aluainua  step  wedge  have  been  constructed  with 
1/8”  steps.  As  the  aluainua  is  about  1/3  as  dense 
as  the  steel,  the  steps  provide  a  24  fold  dynaaic 
range  calibration  in  the  region  of  internet.  These 
wedges  will  be  used  in  every  view  of  exposure.  The 
data  reduction  photodenaltoaeter  will  read  the  wedge 
radiographs  for  each  view  to  provide  a  direct 
density  calibration.  It  is  felt  that  the 
calibration  and  dynaalc  range  error  can  be  kept 
below  10*.  Variation  of  flln  sensitivity  across  a 
single  piece  of  flln  is  felt  to  be  negligible. 

Reading  error  -  a  photodenaltoaeter  has  been 
constructed  and  appears  to  give  repeatable  results 
to  within  about  3*.  Its  accuracy  will  be  further 
evaluated  based  on  our  planned  experlnents. 

Photon  statistics  -  a  radioactive  source  does 
produce  a  tiae  and  space  varying  radiation  field. 

As  we  are  calibrating  every  view,  the  tiae 
canceled  out.  There  is  the  possibility  of  spatial 
varlatlona  but  as  we  are  dealing  with  a  very  large 
nuaber  of  photons,  this  effect  is  negligible. 

Distributed  source  -  the  distributed  source  problea 
is  saall  as  the  sources  in  Bind  are  less  than  1  ca 
in  size  (less  than  the  grid  size).  This  effort  is 
felt  to  be  negligible,  but  will  be  quantified. 


Slicing  error  -  As  discussed  earlier,  the 
radiography  of  any  segnent  that  is  still  connected 
to  the  rest  of  the  body  will  cause  soae  error  due  to 
non-overlap  of  the  varioua  views  at  the  Interface 
between  the  segnent  and  the  rest  of  the  body.  It 
is  felt  that  by  use  of  shallow  angle  views  (but  not 
zero  angle),  this  effect  can  be  kept  below  5*.  This 


6.0  PHOTODENSITOMETER 


A  photodensl toaster  has  been  constructed.  It  allows 
i wasp Ini  a  photocell  over  a  radiograph.  The  position  of  the 
photocell  la  oeasured  using  two  orthogonal  Celesco  lanyard 
transducers  (.008"  accuracy)*  The  transducers  have  ♦  5  volt 
output  for  +  10"  travel.  The  photocell  reads  light 
Intensity  fros  a  reflected  light  source  next  to  the 
photocell.  The  area  of  llluelnatlon  and  sensing  Is 
approxinately  a  1  ce  diaaeter  circle.  The  photocell  output 
is  0  to  8  volts  with  adjustable  sensitivity. 

The  three  signals  generated  (two  coordinates  and  one 
intensity)  are  continuously  read  by  an  IBM- PC  based  digital 
data  acquisition  prograe,  ANPILM.  Bach  1/10  second,  the 
data  Is  sanpled  and  averaged.  The  intenaity  Is  then 
assigned  to  the  geonetric  grid  point  given  by  the  two 
aeasured  coordinates.  The  value  la  also  displayed  on  a  CRT 
screen  grid.  Thus  the  photocell  can  be  swept  by  hand  over 
the  filn  at  randon  without  regard  to  grid  lines  to  "fill  in" 
the  picture.  Moving  slowly  in  significant  areaa  resulta  In 
aultlple  aeasurenenta  and  averaging  of  the  sane  grid  point, 
hence  laprovlng  accuracy.  Once  the  grid  Is  filled  In,  It  is 
written  in  a  standard  data  set  for  future  processing  (called 
an  "R"  data  set). 

ANFILM  also  allows  positioning  of  the  photocell  on  the 
radiographic  wedge  steps  for  calibration.  Thus  the 
photocell  readings  are  directly  and  laaedlately  converted  to 
actual  aass  densities,  cancelling  out  aany  potential  errors. 

AMFILM  also  allows  positioning  on  the  "spot”  points  on 
the  segnent  so  as  to  establish  their  exact  relation  to  the 
aeasureaent  coordinate  systen. 

7.0  SIMULATION  PROGRAMS 

Two  other  prograas,  ANRAY  and  ANSEG,  have  been  designed 
and  are  about  to  be  used.  ANRAY  has  the  purpose  of  solving 
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the  equations  of  section  4.1.2  of  our  proposal.  Thus,  given 
three  "R"  data  sets.  ANRAY  calculates  the  corresponding  aass 
and  inertial  properties.  (Initially  for  Phase  I.  we  are 
assuming  9  orthogonal  views  and  parallel  beaaa.) 

The  second  prograa  ANSEG,  allows  the  user  to 
arbitrarily  define  the  aass  density  of  a  segaent  in  a  3 
diaensional  space  grid.  It  then  calculates  the  aass  and 
inertial  properties  and  the  3  resulting  "R"  data  sets.  This 
prograa  will  be  used  In  the  slaulations  to  evaluate  grid 
size  error,  absorption  variation  (with  atonic  nuaber)  error, 
and  dynaaic  range  requlreaents.  Thus  the  varioue  "R"  sets 
froa  ANSEG  will  be  analyzed  by  ANRAY  to  see  how  well  the 
aass  and  Inertial  properties  can  be  reconstructed  after 
various  error  sources  are  Introduced. 

8.0  PLANNED  EXPERIMENTS 

A  series  of  flla  exposures  are  planned  to  establish 
exposures,  dynaaic  range,  and  accuracies.  The  files  will  be 
exposed  to  both  low  energy  (80  keV)  and  higher  energy  (680 
keV  C137)  radiation  (for  purposes  of  coaparison).  In  all 
exposures,  the  steel  and  alualnua  stepped  wedges  will  be 
used.  Two  phantoas  will  be  eaployed.  One,  for  qualitative 
evaluation,  is  an  approxiaately  8"x8"x8"  plastic  container 
filled  with  water  and  selected  objects  (spotting  ball 
bearings  and  anlaal  bones).  Its  radiograph  will  be  reviewed 
to  deteraine  correct  exposure  and  dynaaic  range.  It  is  aeant 
to  roughly  siaulate  a  priaate  head.  The  second  phantoa  is  a 
6"  homogeneous  plastic  sphere.  Its  radiograph  will  be  used 
to  quantitatively  evaluate  the  ability  of  the  data  reduction 
procedure  to  estlaate  aass  and  inertial  properties. 


